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ABSTRACT 

We report a change in the trajectory of a well studied jet component of the quasar 3C 279. The component 
changes in apparent projected speed and direction, and we find it to be moving with a Lorentz factor 7 > 15 at 
an initial angle of < 1° to the line of sight. The new trajectory of the component has approximately the same 
speed and direction as an earlier superluminal feature, originally seen in the early 1970s. The new direction for 
the component is also much better aligned with larger scale VLBA and VLA structure out to 0.1". We suggest 
that the trajectory change is a collimation event occurring at > 1 kiloparsec (deprojected) along the jet. While the 
change in trajectory on the sky appears to be 26° the intrinsic change is < 1 ° . We estimate the Doppler factor prior 
to the change in direction to be 5 > 28 and after the change to be (5 > 23. Comparison to independent constraints 
on the Doppler factor suggest that the energy in the radiating particles cannot greatly exceed the energy in the 
magnetic field unless the volume filling factor is very much less than one. 

Subject headings: galaxies : active — galaxies: jets — galaxies: kinematics and dynamics — quasars: 
individual: 3C279 



1. INTRODUCTION 

Faster-than-light or "superluminal" motion was originally 
predicted by Rees (1966) and first observed by Whitney et al. 
(1971) and Cohen et al. (1971) in the powerful radio quasar 
3C279 (1253-055; z = 0.536). The illusion of superluminal 
speed results from highly relativistic motion directed toward 
the observer, with the moving material chasing the radiation it 
emits, creating a compression in the time sequence of events as 
seen by the observer. In radio quasars, like 3C 279, a pair of 
highly collimated, relativistic plasma jets stream outward from 
the nucleus or core, and if one of these jets points within a small 
angle, 9, of our line of sight, then pattern motion with some in- 
trinsic speed, j3 = v/c < 1, will appear to move across the sky 
with an apparent speed given by: 
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which can greatly exceed unity. 

The radio jet in 3C279 is a particularly interesting example 
with components moving along multiple position angles and 
apparent speeds. The motion of the original superluminal fea- 
ture, as summarized by Cotton et al. (1979), who analyzed the 
six available epochs spanning the period 1970-1972, was along 
a position angle of-142° ±2° at0.5±0.1 mas/yr. Foramodern 
cosmology of //o = 70kms~^ Mpc"', SIm = 0.3, and J^a = 0.7, 
this corresponds to jSap,, = 16 ±3. The VLA observations of 
de Pater & Perley (1983) showed that knot C at 0.1" had a 
position angle of -145° in a nearly direct line with the super- 
luminal motion detected on length scales fifty times smaller 
The observational picture became less clear during the 1980s 

' National Radio Astronomy Observatory, Charlottesville, VA 22903; dhoman@nrao.edu, mlister@nrao.edu; kkellerm@nrao.edu 
^ Department of Astronomy, MS 105-24, California Institute of Technology, Pasadena, CA 91125; mhc@astro.caltech.edu 

^ Max-Planck-Institut filr Radioastronomie, Auf dem Hilgel 69, D-53121 Bonn, Germany; ros@mpifr-bonn.mpg.de, azensus@mpifr-bonn.mpg.de, mkadler@mpifr- 
bonn.mpg.de 

Netherlands Foundation for Research in Astronomy, Postbus 2, NL-7990 AA Dwingeloo, Netherlands; rvermeulen@astron.nl 

' The VLBA is operated by the National Radio Astronomy Observatory which is a facility of the National Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 

1 



when Unwin et al. (1989) tracked another superluminal compo- 
nent from 3C 279 which they called "C3". Component C3 was 
found to move along a different position angle of -134° with a 
much smaller speed of only 0. 12 ± 0.02 mas/yr corresponding 
to/3flpp=4±l. 

At the very end of 1984, Unwin et al. (1989) note the emer- 
gence of a new component they call C4. Component C4 is an 
extremely strong and compact feature which has been tracked 
up to the present by a number of programs. The most detailed 
and complete results to date were presented by Wehrle et al. 
(2001). They followed C4 from 1991 to the end of 1997 with 
22 GHz observations and from 1995 to 1997 at 43 GHz. They 
found the component to be moving along a position angle of 
-114° with a speed of 0.26 ±0.01 mas/yr which corresponds 
to 8.2 ± 0.3 times the speed of light in our choice of cosmol- 
ogy. They extrapolate an origin epoch for the component of 
1984.7 ±0.3, consistent with the first observation of C4 by Un- 
win et al. (1989). Wehrle et al. (2001) discuss a dip in speed 
of C4 around 1994 which coincides with small changes (a few 
degrees) of the trajectory on the sky. Whether this small change 
in the motion is related to the large changes reported here is un- 
clear; however, the component had resumed its previous motion 
and trajectory by the end of 1995. 

Here we present results from the Very Long Base- 
line Array^ (VLBA) 2cm Survey (Kellermann et al. 1998, 
http://www.cv.nrao.edu/2cmsurvey) showing that after 1998, 
the trajectory of C4 changes in both apparent speed and direc- 
tion. The new trajectory for C4 is on an essentially parallel 
track to the 1970s superluminal component with essentially the 
same speed as originally observed for that component. 
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Fig. 1 .— a A2 cm, VLBA image of 3C 279 in Marcli 1997. Contours begin 
at 10 mjy/beam and increase in powers of two. The cross in the NW comer of 
the map depicts the FWHM dimensions of the restoring beam. Overlayed on 
this image in red are component positions for C4. Solid triangles are 1 1 GHz 
data from Carrara et al. (1993), open squares are 22/43 GHz data from Wehrle 
et al. (2001) (pre 1995.5), and solid circles are measured from 2 cm VLBA 
data taken from the period 1995.5 to 2002.9. Blue open star's and green open 
crosses indicate the approximate locations measured for the 1970s superlumi- 
nal component and component C3 respectively (distances and position angle 
taken from Cotton et al. (1979) and Unwin et al. (1989)). Red, green, and blue 
dashed lines mark the —1 14°, -134°, and —142° position angles respectively. 



2. OBSERVATIONS AND RESULTS 

Our observations of 3C 279 from the VLBA 2cm Survey con- 
sist of 12 epochs spanning the interval from 1995.5 to 2002.9 
(for an animation see our project website referenced above). 
We have also incorporated 2cm VLBA results from the Bran- 
deis University parsec-scale jet monitoring project (Homan et 
al. 2001, 2002, Ojha et al. in prep.) which adds six closely 
spaced epochs during 1996 and a seventh epoch at the end of 
1997. We measured component positions both in the image 
plane, by fitting point sources at the peak locations of the core 
and C4, and in the (u, v)-plane by fitting 2 or 3 point compo- 
nents to the core region and a single elliptical Gaussian to C4. 
We have averaged the values obtained by these two approaches, 
and have plotted the mean position in figures 1 and 2. The core 
is assumed to be stationary, and the positions of C4 are taken 
relative to the core. 

Figure 1 shows an image of 3C279 from our March 1997 
observations. Superimposed on this image is the position of 
component C4 from each of the 19 epochs we compiled. We 
have also plotted 22/43 GHz data from Wehrle et al. (2001) up 
to 1995.5 and 11 GHz data from CaiTara et al. (1993) to illus- 
trate the trajectory of C4 from 1988 to the present. Approxi- 
mate positions for the 1970s superluminal component and the 
component C3 are included on the figure for comparison. 

Proper motions were computed from the 2cm (x,y) compo- 
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Fig. 2. — Component locations and fitted proper motions of C4 before and 
after the bend. Panel (a) plots the component locations and fitted motions as 
seen on the sky. Panel (b) plots the component distances and fitted motions 
versus time. The distances in panel (b) are projected distances along the po- 
sition angles con'esponding to the directions of motion before and after the 
bend. 



nent positions in both the image and (m, v) planes independently. 
The motions were computed over two different intervals, one 
prior to the change in trajectory and one after We were unable 
to decide whether epoch 1998.2 should be classified with the 
motion prior to the change or after the change, but we found 
that it made little difference which choice we made. The aver- 
age motions we found across all of these approaches were the 
following: priorto ^1998.2: ^ = 0.25±0.02 mas/yr (7.9 ±0.6 
c) along a position angle of -114° ± 1°, and after ^1998.2: 
fi = 0.40 ± 0.01 mas/yr (12.7 ± 0.3 c) along PA= -140° ±1°. 
These proper motions are plotted against the data in figure 2. 
We estimate our typical positional uncertainty for the centroid 
of C4 to be 0.03 mas in both dimensions. 

Figure 3 depicts the change in component flux density (panel 
(a)) and size along the -114° position angle (panel (b)) as a 
function of epoch. Prior to the sharp change in motion some- 
time in 1998, the component increases in flux density by 50% 
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Fig . 3 . — Plots of component C4's flux density (panel (a)) and FWHM angu- 
lar size projected along the —1 14° position angle (panel (b)). These properties 
are taken from our Gaussian fit to the (u,v)-plane data. The square symbols 
represent a poorly calibrated epoch with flux density ~ 10% low. Error bars 
on the flux density points are estimated at 5% of the integrated flux density in 
each epoch (e.g. Homan et al. (2002)). 

and decreases in size along the direction of motion by a factor 
of two. These events seem to come before the change in mo- 
tion, suggesting the component is already undergoing changes, 
likely caused by interactions with its environment, while it is 
still propagating along the original direction of motion. Zavala 
& Taylor (2001) report a change in the polarization position 
angle of C4 which occurred sometime between 1998.6 (Tay- 
lor 2000) and 2000.1. The polarization angle change reflects a 
change in magnetic field order, as Faraday rotation is small in 
this feature (Zavala & Taylor 2001), and provides further evi- 
dence that the change in trajectory is tied to internal changes in 
component structure. 

3. ANALYSIS 

Here we make a simple "billiard-ball" analysis of the com- 
ponent's trajectory. We note that more complicated models 
are possible, perhaps involving very different pattern and flow 
speeds, and such models may allow a wider range of aUowed 7 
and e. 

The apparent speed of this component has increased by more 
than 50% and the trajectory on the sky has changed by 26°. 
According to equation 1 , the apparent speed can increase either 
due to an increase in the intrinsic speed, /3, or by a change in an- 
gle to the Une of sight. We observe a large change in projected 
angle on the sky, therefore it is reasonable to assume that the 
angle to the line of sight has also changed. We cannot rule out 
an increase in beta; however, the simplest assumption consis- 
tent with our observations is that the speed of the component is 
constant along a bent path, and we make that assumption here. 

Figure 4 plots the Doppler factor, 6 = 1 /7(1 -/3cos^) where 
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Fig. 4. — Plot of Doppler factor versus apparent speed. Dashed contours in- 
dicate lines of constant bulk Lorentz factor, and solid contours indicate lines of 
constant angle to the Une of sight. Vertical lines with gray envelopes indicate 
the component speed before (7.9 ± 0.6 c) and after (12.7 ± 0.3 c) the bend. 

7=1/ 1/1 versus the apparent motion, papp- Lines of con- 
stant speed (Lorentz factor, 7) and lines of constant angle are 
shown on the figure. 7 = 10 is clearly too small and cannot ac- 
commodate speeds greater than 10c. The lines for 7 = 15 or 
7 = 20 are more reasonable, and the change in apparent speed 
from 8c to 1 3c can be accommodated by either a smaU increase, 
~ 1 °, or a larger decrease, 6°, in the angle of the motion with 
respect to the line of sight. 

To decide between the two possibilities of increasing or 
decreasing 6, we must look at the corresponding change in 
Doppler factor which will scale the flux density of the jet com- 
ponent like S oc S^-^, appropriate for a simple discrete compo- 
nent with a spectral index of 0.5. For a decreasing 9 (i.e. a 
trajectory bent towards the line of sight), the Doppler factor 
more than doubles between 8c and 13c along the 7 = 15 curve. 
This would predict an increase in flux density of the component 
of more than an order of magnitude; however, from figure 3a, 
it is clear that the flux density of the component changes by no 
more than a factor of two, and after 1998 the flux density only 
decreases. Note that the observed compression in component 
size prior to the bend would be likely only to further increase 
the flux density, and therefore with some confidence we reject 
a decreasing angle to the line of sight. For increasing 9 (i.e. 
a trajectory bent away from the Une of sight) the situation is 
much better as the Doppler factor decreases by approximately 
20% along the 7 = 15 curve. This corresponds to a decrease in 
component flux density of a factor of two, somewhat more than 
we actually observe but not unreasonably so. 

Considering the complex nature of the light curve which is 
influenced not only by Doppler beaming but also by the ob- 
served compression of the component prior to the bend and 
expansion after the bend, we don't place strict limits here on 
how the Doppler factor has decreased after the bend. While 
the 7=15 curve seems to provide a reasonable lower bound on 
C4's Lorentz factor, we note that the Lorentz factor could be 
much larger and still be consistent with our observations. 

From this analysis assuming fixed component speed, we es- 
timate 7 ^ 15 with an initial angle to the line of sight of 6 <l° 
increasing at the bend to become 6 <2°. The bend in the plane 
of the sky appears to be 26°; however, deprojected, this bend 
is only ~ 26° x sin^ ~ 0.5° - 1°, consistent with the required 
bend along the line of sight. 
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4. DISCUSSION 

4.1. The Nature of Component C4 

In the above analysis, we have taken a very simple "billiard- 
ball" model of the motion of component C4, where C4 under- 
goes a sudden change in trajectory sometime during 1998. With 
a major axis size (FWHM) of < 10 ly (0.5 mas) and a Doppler 
factor > 20, it is not unreasonable that C4 could have changed 
its trajectory in less than a year in our frame, and while this 
simple model fits our data quite well, a more detailed analysis 
of the trajectory is desirable. Unfortunately we are limited in 
this by a lack of good sampUng in the period 1998-2000. An 
additional problem for more detailed analysis of the trajectory 
is that C4 is an extended region, and we are only tracking its 
centroid. In reality, C4 may be a complex shocked region with 
a strong leading edge and trailing shocks, and small changes in 
its centroid may not reflect true kinematical changes. In this 
light, we believe our simple approach is a conservative one that 
leads to a consistent picture of the jet kinematics. 

4.2. Doppler Factor and Equipartition 

From the angles and Lorentz factors estimated above, we can 
derive approximate hmits on the Doppler beaming factor of this 
jet feature before and after the bend. We find ^ > 28 prior to 
the change in trajectory and 5 > 23 after the change. These val- 
ues are consistent with Finer et al. (2003), who use an inverse- 
Compton calculation to estimate an average lower limit for the 
Doppler factor of C4, ^ > 21, prior to the change in direction. 

It is interesting to compare our values to the completely inde- 
pendent estimate of the Doppler factor made by Homan & War- 
die (2000) assuming equipartition between the radiating parti- 
cles and magnetic field. They observed a spectral turnover at 6 
GHz for C4 in epoch 1997.94, just before the change in trajec- 
tory described here. They used the spectrum and resolved angu- 
lar size of C4 to calculate the Doppler factor: 5 = IS^^^v/^)^^'' 
where rj is the ratio of energy in the magnetic field to energy 
in the radiating particles: ry = Us/Up, and where T < 1 is the 
fraction of the volume filled with radiating plasma^. Agreement 
between these two independent determinations of the Doppler 
factor of component C4 requires ry/T > 1, implying that the en- 
ergy in the radiating particles cannot greatly exceed the energy 
in the magnetic field unless the volume fiUing factor is very 
much less than unity. 



4.3. What Caused the Bend? 

We argue above that C4 was interacting with its environment 
prior to the change in trajectory, resulting in decreased angular 
size along the -1 14° position angle and increased flux density. 
It seems likely, therefore, that C4 has been deflected onto its 
new trajectory. A random collision would be unlikely to de- 
flect C4 in just the right fashion to have essentially the same 
speed and direction on the sky as the superluminal component 
detected in 1971. It is important to stress that while the trajec- 
tory of C4 is parallel to that of the older component, C4 is not 
following the actual path of that component (see figure 1). The 
new trajectory for C4 not only has a very similar speed and di- 
rection as the 1970s superluminal feature, but it is also in the 
direction of the larger scale structure as seen in the 1.6 GHz 
VLBA maps made by Finer et al. (2000). 

We suggest that the change in trajectory is a collimation event 
resulting from the interaction of C4 with the boundary between 
the jet outflow and the interstellar medium. The exact nature of 
this boundary and interaction is unclear; however, it is interest- 
ing to note that (a) the apparent collimation event occurred at a 
deprojected distance from the nucleus of > 1 kpc, and (b) the 
jet is still highly relativistic, with 7 > 15, on these scales. 

Abraham & Carrara (1998) discuss a precessing jet model for 
3C 279 to explain the different apparent speeds and position an- 
gles of the VLBl components appearing prior to 1990. Wehrle 
et al. (2001) confirm the presence of different apparent speeds 
and position angles in the jet, even in newer components origi- 
nating after C4, although their measurements do not agree with 
the extrapolated predictions of the Abraham and Carrara model. 
The concordance of the new speed for C4 with that of the 1970s 
component adds weight to the idea that jet components may all 
have intrinsically similar Lorentz factors but shghtly different 
angles of ejection. If some sort of precession model is used 
to interpret the data, Wehrle et al. comment that it is difficult 
to understand how the range of ejection angles from -114° to 
-142° can be consistent with the narrow range of angles ob- 
served in the larger scale VLBA and VLA structure. Ferhaps 
we have seen the answer here in the form of jet colhmation that 
occurs at ^ 1 kpc or more. 

We thank R. Ojha, D. Roberts, and J. Wardle for sharing 
their monitoring data. J. A. Z. was supported for this research 
through a Max-Flanck Research Award. This research has 
made use of data from the University of Michigan Radio As- 
tronomy Observatory, which is supported by the National Sci- 
ence Foundation and by funds from the University of Michigan. 
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